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Abstract 

Many  service  systems  confront  tasks  of  uncertain  identity,  with  limited 
time  for  service;  that  is,  the  uncertain  tasks  have  deadlines,  or  are 
behaviorally  impatient.  Examples  occur  in  medical  care  (especially 
emergencies),  telephone  help  systems,  and  in  military  operations.  This 
paper  presents  modifications  of  the  M/G/l  system  to  illustrate  the  impact 
of  the  above  features.  Imperfect  task  classification  is  modeled,  as  is 
imperfect  service  and  error-afflicted  assessment:  tasks  can  be  processed, 
and  reprocessed,  either  correctly  or  incorrectly  depending  upon 
classification,  performance,  and  performance  assessment  skills.  The 
impact  of  exponential  deadlines,  either  behavioral  or  server-controlled,  is 
represented  using  both  a  modification  of  the  Takafs-Benes  integro- 
differential  equation,  and  a  simple  and  accurate  fixed-point  approximation. 

1.  Background 

Many  service  situations  are  characterized  by  currently  unmodeled  uncertainties  at 
least  as  influential  as  those  identified  with  the  usual  stochastic  arrival  and  service 
processes.  Additional  and  important  uncertainty  sources  or  components  include  the  true 
identity  or  nature  of  the  task,  hence  the  “optimal”,  or  at  least  satisfactory,  mode  of  its 
service  (in  the  light  of  resources  available  and  the  other  tasks  on  hand  or  anticipated). 
Further  complications  occur  if  the  task  has  a  deadline:  is  impatient  or  perishable,  i.e.  has 
an  unknown  life  or  time  available  for  useful  service:  e.g.  telephone  callers  placed  on  hold 


(prone  to  hang  up  or  abandon  an  attempt),  medical  emergency  patients  (who  may  die), 
mobile  military  targets  (that  tend  to  move,  or  possibly  fire  first),  and  many  others. 

These  questions  naturally  arise:  once  a  service  is  performed,  is  it  “complete”,  or 
should  additional  service  work,  possibly  of  a  different  kind,  be  performed?  And  what 
fraction  of  time-sensitive  perishable  or  deadline-afflicted  tasks  actually  finish  service? 
Such  issues  arise  because  the  initial  classification  of  tasks  is  realistically  uncertain,  as  is 
the  degree  of  task  accomplishment.  Post-service  assessment  is,  also  realistically,  error- 
prone,  which  can  sometimes  lead  to  premature  release  of  an  incompletely  served  task, 
with  costly  legacy,  or  to  the  needless  expenditure  of  extra  time  and  resources  to 
“complete”  a  task  already  completed.  We  call  the  above  generic  uncertain  service 
situations,  and  study  several  of  their  features.  Important  practical  issues  are  to  discover 
system  performance  sensitivities,  such  as  which  additional  uncertainty-reduction 
capabilities  are  likely  to  be  most  cost-effective. 

As  emphasized  above,  uncertain  service  situations  arise  in  many  contexts.  Examples 
are  medical  emergencies,  e.g.  those  that  arise  from  terrorist  attacks  or  natural  hazards, 
such  as  hurricanes,  floods,  or  earthquakes;  forest-fire  fighting;  engineering  and 
operational  problems  with  new  computer  software  and  hardware;  and  in  military 
operations,  e.g.  on  the  battlefield,  where  “service”  means  actual  or  potential  destruction 
or  deterrence  of  opponent  assets  or  troops,  but  may  also  involve  inadvertent  server 
deception,  fratricide,  collateral  damage,  and  injury  to  non-combatants.  Avoidance  of 
these  latter  limits  service  options. 

In  this  report  the  deadlines,  or  impatience,  that  characterize  our  tasks  are  described 
probabilistically:  the  server  only  knows  the  distribution  of  deadline  elapse  and  task 
disappearance,  not  individual  durations.  There  is  a  large  literature  on  task  scheduling  in 
the  face  of  hard,  known,  deadlines;  cf.  Liu  and  Layland  (1973),  Jiang,  Lewis,  and  Colin 
(1996),  and  others.  More  recently,  stochastic  scheduling  of  queued  tasks  with  deadlines 
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that  are  subjected  to  particular  queue  disciplines  are  studied;  cf.  Lehoczky  (1996,  1997a, 
1997b),  and  Doytchinov,  Lehoczky  and  Shreve  (1998).  The  latter  problems  come  from 
situations  for  which  it  is  natural  to  know  the  deadline;  our  examples  are  otherwise. 
Imperfect  service  (repair)  has  also  been  studied  by  Brown  and  Proschan  (1983),  but  not  in 
a  congested  setting.  Quite  possibly  there  is  related  work  by  others  that  is  unknown  to  us. 
A  specialized  and  more  intricate  deterministic  version  of  our  problem  is  described  in 
Gaver  and  Jacobs  (1999).  Software  illustrating  that  model  is  available  from  the  authors. 

2.  Analytical  Setting 

The  purpose  of  this  paper  is  to  explore  and  expose  operating  characteristics  of  a 
simplified  totally  stochastic  version  of  the  above  setup.  The  mathematical-probability 
theory  of  a  single-server  queue,  cf.  Cox  and  Smith  (1961),  Kleinrock  (1976),  and  many 
others,  is  adapted  to  study  an  uncertain  service  situation  in  which  deterioration  of  task 
value  also  occurs  (deadlines  are  missed,  patients  die,  military  targets  move). 

Our  basic  model  utilizes  a  first-in,  first-out  (FIFO)  basic  queue  discipline  for  several 
reasons.  Simplicity  is  first  among  these:  understanding  how  to  establish  priorities  in  the 
face  of  the  other  uncertainties  faced  is  difficult,  and  is  postponed.  Real-time  maintenance 
of  priority  or  other  control  is  also  time-consuming,  so  the  time  cost  of  control  should  be 
included;  see  the  model  of  Appendix  C  for  a  start.  Alternatively,  one  can  view  the  present 
model  as  being  that  for  a  service  system  that  experiences  triaged  traffic  from  an  initial 
screener.  Disciplines  other  than  FIFO  are  ignored  at  this  level.  We  do  plan  to  attack  the 
uncertain  service  problem  in  a  more  general  control  environment  in  future. 

3.  Model  Formulation 

A  task  of  type  j  is  classified  as  being  of  type  k  with  probability  cjk;  to  be  interesting, 
Cjj  <  1,  meaning  that  task  misclassification  may  occur,  perhaps  with  appreciable 
probability.  A  model  input  is  the  classification  matrix  c.  The  (possibly  misclassified)  task 
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is  then  prosecuted  with  probability  of  success  mjk;  this  probability  depends  on  treatment 
choice;  let  mjk  =  1  -  mjk  be  the  probability  of  failure.  A  model  input  is  success  probability 

matrix,  m.  Assessment  (service  outcome  inspection)  is  next  conducted;  if  the  task  has 
been  completed/carried  out  successfully  this  fact  is  ascertained  correctly  (verified)  with 
probability  bjk;  if  the  task  has  been  completed  successfully  it  is  incorrectly  reported  as  not 
complete,  hence  is  a  candidate  for  needless  repetition  with  probability  1  -  bjk  =  bJk.  If  the 

task  is  not  completed  successfully  it  is  correctly  reported  as  non-successful  with 
probability  b]k,  and  incorrectly  reported  as  successful  with  probability  b]k  =  1  -b*jk.  All  of 

these  parameters  can  be  made  functions  of  other  variables:  for  example,  patient  condition 
such  as  age  or  auxiliary  measures  and  symptoms  in  medical  diagnosis,  or  such  as  range, 
atmospheric  conditions,  terrain,  deception  tactics  in  military  command  and  control.  We 
treat  them  as  constants  for  the  present.  Thus  further  model  inputs  are  the  post-service 
assessment  matrices  b  and  b*.  The  actual  time  spent  servicing  the  task,  which  depends 
upon  the  classified  state,  is  subject  to  the  decision  maker’s  influence  (e.g.  it  may  be 
truncated);  it  too  is  a  decision  variable. 

In  general,  there  may  be  a  number  of  attempts  made  to  complete  a  task,  and  to 
confirm  that  completeness.  Depending  on  the  classification  (c),  the  probability  of  success 
and  task  completion  (m),  and  assessment  (b,  b*)  skills  of  the  service  facility,  the  system 
will  either  provide  good  and  timely  service,  or  not  five  up  to  its  promise.  In  various  cases, 
the  reason  for  performance  degradation  may  well  depend  dominantly  upon  timely 
information  available  concerning  the  task  and  its  accomplishment,  and  less  on  the  actual 
or  true  probability  of  task  prosecution  success.  In  all  real  situations  the  required 
information  and  capabilities  are  only  available  at  a  cost.  It  is  the  purpose  of  this  study  to 
illustrate  the  nature  of  the  various  possible  cost  tradeoffs.  The  models  proposed  are  a 
beginning,  but  ultimately  a  means  to  that  end. 
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4.  Stochastic  Model 

Arrivals  appear  at  a  service  facility  according  to  a  Poisson  (A)  process.  The 
probability  that  an  arrival  is  of  type  j  is  pj  (j=  l,  2,  independently  from  task  to  task. 
Let  Sk  denote  the  processing  time  initially  allotted  to  a  task  classified  as  of  type  k,  an 
assigned  service  time.  Note  that  this  is  not  necessarily  (or  even  frequently,  in  the  present 
context)  the  time  to  successfully  service  an  item  of  true  type  j,  especially  one  different 
from  type  j.  We  use  Sk  to  represent  the  time  to  carry  out  a  particular  process  that  it  has 
been  selected  to  apply,  distinguishing  this  from  the  probability  of  process  success,  mJk,  or 
mjk(Sk)  if  desired.  As  suggested,  Sk  may  have  a  decision  component,  i.e.  be  subject  to  a 
decision  maker’s  choice. 

4.1  Individual  Server  Occupancy  Times  (ISOT) 

Suppose  a  task  of  type  j  presents  itself  to  a  servicing  facility.  In  what  follows  we 
characterize  its  continuous  occupancy  of  the  servicing  facility,  e.g.  a  diagnostic  and 
treatment  sojourn  with  a  medical  facility,  or  as  the  current  target  of  a  generic  shooter  in  a 
military  context.  The  task  may  actually  complete  long  before  the  generalized  server 
recognizes  that  fact;  on  the  other  hand,  the  generalized  server  (server  plus  reassessment 
asset)  may  act,  and  prematurely  and  incorrectly  decide  that  the  task  is  complete.  Some 
“completed”  tasks  are  thus  released  in  misdiagnosed  and  dangerous  condition,  either  to 
themselves  (medicine)  or  others  (military).  Our  models  allow  understanding  of  system 
tradeoffs  that  control  the  probability  of  such  happenings. 

4.2  Random  Reclassification  after  Each  Assessment 

Suppose  the  system  is  arranged  so  that  reclassification  occurs  independently  and 
“with  replacement”  immediately  after  each  assessment  that  declares  an  unsuccessful 
service  attempt.  If  the  assessment  declares  success  a  new  task  begins.  This  is  just  one 
simple  option ;  see  Appendix  C,  which  proposes  a  decision  rule  that  may  reclassify  if  the 
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task  service  is  perceived  to  be  incomplete;  more  sophisticated  options  are  available  at  a 
price  in  time.  Here  is  the  corresponding  model  for  the  individual  server  occupancy  time 
(ISOT)  of  a  task  of  true/actual  type  j  under  random  reclassification: 
for j,k=  1,2,  ...,J 


Cj  = 


sk+q 

Sk  +  Kj 


with  probability  cjk  [m jk  (Sk  )bjk  +  mjk  (Sk  )bjk  ]; 
with  probability  cjkmjk(Sk)bjk ; 
with  probability  cjkmjk(Sk)bjk; 


(4.1) 


Kj  represents  the  random  time  until  a  completed  task  is  so  identified,  the  recognition  time : 
'Sk  with  probability  cjk  •  1  ■  bjk 

Kj=\  _  (4.2) 

Sk  +  K}  with  probability  cjk  •  1  •  bJk 

Notice  that  the  task  accomplishment  probability  is  allowed  to  depend  explicitly  on  the 
allocated  service  time,  Sk,  and  that  the  ISOT  can  terminate  with  unrecognized  incomplete 
task  service,  i.e.  the  task  may  be  terminated  although  incompletely  served.  The  random 
variables  C'  and  K)  above  are  independent  stochastic  replicas  of  Cj  and  Kj. 


4.3  Expectations 

Taking  conditional  expectations  we  obtain  these  expressions  for  mean  ISOT: 

£[C,]=y^£[5,]  +  y  (St  )]t^] + £  <*£[»,(*  )JAKA  <43) 

k  k  k 

where 

E[Kj]  =  ^cjkE[Sk]  +  ^cjkbjkE[Kj],  (4.4) 

k  k 


which  leads  to  the  formula 


l  +  y \-'YJCjUb,k 


1  ^  Cjkkfytljk  i$k  'ityjk 


(4.5) 
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An  expression  for  the  second  moment  appears  in  Appendix  A.  It  is  very  clear  from  (4.5) 
that  degradation  of  performance  can  be  associated  with  misclassification  and  faulty 
assessment. 

4.4  Task  Queue 

Tasks  appear  at  the  service  facility  at  Poisson  rate  X.  Any  task  is,  independently,  of 
type  j  with  probability pj  (j  =1,2, ...,  J),  but  the  type  is  only  known  with  uncertainty,  cjk. 
If  tasks  are  treated  according  to  the  first-come,  first-served  discipline  by  a  single  server 
then  the  system  is  M/G/l  with  effective  service  times 

e\.c\=y.pAcA  (4-«) 

j 

hence,  traffic  intensity 

p=XE[C ].  (4.7) 

One  measure  of  the  system  congestion  is  then  the  long-run  expected  number  of  enqueued 
tasks  (if  p<  1): 

(4.3) 

This  measure  does  not  reflect  the  number  of  tasks  that  are  terminated  before  service  is 
complete.  If  a  goal  is  to  minimize  weighted  expected  waiting  time  then  prioritization  in 
accordance  with  the  index  wJE[Cj\  is  optimal.  Here  Wj  is  the  desirability  weight 
associated  with  completing  class  j;  task  groups  are  served  in  order  of  increasing  index. 
This,  however,  takes  no  account  of  the  influence  of  deadlines.  When  deadlines  are  an 
important  feature  of  the  problem  then  different  measures  of  system  effectiveness  are 
needed. 


7 


5.  Successful  Accomplishment  Probability  for  Tasks  with  Deadlines 

In  a  variety  of  contexts  the  value  of  task  service,  or  the  probability  of  successful 
completion,  decreases  with  the  delay  experienced.  This  is  often  true  of  medical  diagnosis 
and  treatment,  particularly  in  emergencies,  and  in  military  targeting  (the  object  targeted 
may  move).  Using  adaptations  of  the  M/G/l  queuing  model,  cf.  Cox  and  Smith  (1961), 
we  study  tradeoffs  among  service  capabilities  in  such  situations.  Queuing  models  in 
which  an  arriving  customer  is  lost  when  it  waits  more  than  a  fixed  time  in  queue  have 
been  studied  by  Boots  and  Tijms  (1999),  and  Whitt  (1999).  Note  that  server-imposed 
deadlines  are  a  control  device  that  may  improve  certain  measures  of  system  performance. 
These  are  in  effect  a  refusal  to  provide  service;  balking  has  traditionally  been  a  task- 
initiated  refusal;  see  Whitt  (1999).  Our  models  address  refusals  in  general. 


5.1  The  M/G/l  Service  System  with  Ignored  Exponential  Deadlines 

Express  the  delay  sensitivity  or  deadline  for  tasks  of  type  j  arriving  at  the  M/G/l 
system  above  as  an  exponential  random  variable  with  rate  A  well-known  queuing 
theory  result  is  that  the  long-run  waiting  time,  W,  in  an  M/G/l  system  has  Laplace- 
Stieltjes  transform 


E[ e~sW]  = 


1  -P 


1  -E 


sE\ 


c] 


_  1  -P 

i-M5) 


(5.1) 


where  the  role  of  service  time  is  played  here  by  the  ISOT,  C.  For  an  arriving  task  of  type  j 
her  probability  of  surviving  the  wait  in  queue  without  deadline  elapse  is  e~e‘w , 
conditional  on  W,  so,  unconditionally,  the  probability  of  initial  wait-survival  is 

1-/7 


E[e~^w]  = 


l- A 


\-E[e'6jCY 


(5.2) 
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See  Appendices  B  and  C  for  expressions  for  E[e~sC].  In  the  present  model  all  tasks  are 
served  to  completion,  regardless  of  deadline  elapse.  This  may  be  more  reasonable  in 
some  situations  than  others;  it  is  changed  in  a  subsequent  model. 

By  the  memoryless/Markov  exponential  property  the  type  j  task  survives  a  subsequent 

completion  time,  duration  Dp  that  terminates  with  successful  service  with  probability 
E^e~djDj  J ;  see  Appendix  B  for  the  transform  of  the  improper/dishonest  random  variable 

It  follows  that  the  long-run  marginal  probability  that  a  random  task  completes  service 
satisfactorily  is 

.P(paraineters)  =  •  E^e~8jDj  ] 

(5.3) 

e-™} 


The  above  analytical  expression  may  be  evaluated  numerically,  and  explored  for 


parameter  dependencies.  The  results  of  such  investigations  appear  later. 


Note  that  the  above  model  does  not  assume  the  capability  of  detecting  “dead”  or 
deadline-elapsed  tasks  in  the  queue  or  upon  entering  service.  Under  many  conditions  such 
could  be  refused  or  purged,  thus  increasing  the  chance  of  successful  service  for  others. 
The  above  model  results  thus  tend  to  be  pessimistic  or  conservative  from  the  server 
perspective,  but  not  necessarily  unrealistically  so:  additional  capabilities  may  be  needed, 
but  unavailable,  to  monitor  enqueued  tasks  for  real-time  viability.  The  next  model 
addresses  such  capability  by  the  server. 
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5.2  The  M/G/l  System  with  Deadline-Sensitive  Delay 

Consider  the  arrival  of  tasks  with  an  exponential  (0)  deadline,  and  suppose  that  when 
the  task  reaches  the  end  of  the  queue  it  (the  probability)  can  be  determined  that  the 
deadline  will  not  elapse  before  reaching  the  server,  given  the  virtual  waiting  time,  W(t). 
With  that  probability,  the  task  is  accepted  into  the  queue.  We  first  propose  the  following 
heuristic  analyses,  but  follow  up  with  a  more  formal  treatment  in  Appendix  D. 
Approximation  I 

Given  the  waiting  time  encountered  on  generic  arrival,  W,  the  effective  service  (ISOT) 
time  is 
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This  expression  asserts  that  the  probability  of  successful  (deadline  unviolated)  task  arrival 
for  service,  jA.9)-,  is  that  of  an  M/G/l  system  whose  arrivals  are  filtered  by  the  same 
probability  in  the  long  run.  In  effect,  each  arrival  flips  a  biased  coin  with  success 
probability  E\e~ew\  to  be  permitted  to  join  the  queue.  The  result  differs  somewhat  from 
the  solution  of  the  forward  Kolmogorov  (Taka?s-Benes)  equation  for  the  same  assumed 
arrival-queue  interaction;  see  Appendix  D. 

The  expression  (5.6)  is  a  quadratic  in  the  desired  probability,  the  solution  of  which  is 

- r=T - '  (5-7) 

where 


m= 


1  -£[e-*] 
0E[C\ 


(5.8) 


the  transform  of  the  service/completion  time  tail  or  survivor  distribution.  Of  course  the 
probability  of  successful  transit  to  the  server  is  unity  when  0  (no  degradation,  or 
infinite  deadline),  regardless  of  the  value  of  p<  1;  if  oo  then,  since  deadlines  are  now 
stringent,  the  only  hope  of  initiating  service  is  to  arrive  when  there  is  no  server  activity, 
i.e.  with  probability  1/(1  +  p),  irrespective  of  (positive)  /7-value.  Likewise,  there  are  no 
restrictions  on  p  in  (5.7):  a  long  queue  generates  many  rejections,  and  thus  does  not 
remain  long,  or  grow  indefinitely.  Numerically,  the  above  simple  expressions,  (5.7)  and 
(5.16),  supply  a  lower  bound  that  has  been  shown  numerically  to  be  a  good 
approximation  to  the  exact  solution  of  a  refusal  model  proposed  in  Appendix  D.  Note  that 
the  present  approximation  gives  for  the  transform  of  virtual  waiting  time  of  non-refused 
tasks,  W,  the  formula 


1  -Aip(0)E[C 

1  -Aip[0)E[C} 

l-E 

V*]' 

& 

M  | 

(5.9) 
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where  y/^0)  is  given  by  (5.7).  This  can  be  modified  to  represent  refusal  during  service, 
and  to  provide  approximations  to  the  long-run  mean  waiting  time.  No  numerical 
discussion  or  comparisons  are  available  at  present,  although  comparison  to  results  from 
(D.12)  and  (D.13)  is  of  interest.  These  models  can  be  compared  to  those  of  Whitt  (1999). 

Approximation  II 

A  refined  version  of  the  above  accounts  for  the  different  experience  of  a  new  task  that 
arrives  to  find  the  server  busy  (W>  0),  as  contrasted  to  one  that  arrives  to  find  it  idle 
(W=0).  Put 


y/+{&)  =  E[e*y\W>0] 


(5.10) 


the  marginal  long-run  rate  of  task  acceptance  given  that  the  server  is  busy.  From  (5.1) 

E[e~sW\W>0]  = 


i -pKs) 

Approximate  as  follows  (5  =  6=  1/mean  deadline) 

¥A^-P¥A^x_p^me)’  (512) 

this  asserts  that  the  probability  that  an  arriving  task  that  encounters  a  busy  period  and 
reaches  the  service  stage  is  that  of  an  M/G/l  system  whose  busy-period  arrivals  are 
filtered  by  the  same  probability.  In  other  words,  an  auxiliary  randomization  (biased  coin 
flip)  adjusts  for  the  imposition  of  the  deadline,  as  before  in  Approximation  I,  but  in  a 
somewhat  more  refined  manner.  The  solution  of  (5.1 1)  is 


1  -P 

1  ~MS). 
(1 


P 


(5.11) 


vAd)= 


_ 2S{9) _ 

(l + pS[0))  +  -^(l +pS(0))  -ApS[9) 


(5.13) 


For  such  a  ^-filtered  system  the  expected  duration  of  a  busy  period,  E[E],  satisfies 
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(5.14) 


y/^9)  =  P{W  =  0} +(1  -  P{W  =  0  })y+{0) 
1 

l+4l 


(5.16) 


which  differs  from  (5.7)  owing  to  the  more  refined  conditioning  imposed. 

It  will  be  seen  that  Approximation  II  improves  on  Approximation  I  in  all  cases 
considered. 


6.  Numerical  Exploration 

In  this  section  we  display  graphs  of  the  effects  of  the  various  capability  parameters  on 
long-run  probability  of  successful  service  completion. 

For  reference,  the  following  parameters  are  varied 

•  A;  the  Poisson  arrival  rate  of  tasks; 

•  cjk;  (j,k  e[l,  2,.../|):  probability  of  initially  classifying  a  task  of  type  j  as  being 
type  k; 

•  mjh  or  generally  mjk(S *):  probability  of  success  of  a  service  of  task  type  j  when 
prosecuted/served  as  type  k;  mjk  =  1  -  mJk  is  the  probability  of  an  unsuccessful 

outcome; 

•  bjk:  probability  of  successful/correct  assessment  as  successful  treatment  of  task  of 
type  j,  given  it  has  been  successfully  prosecuted/served  as  type  k;  1  -  bJk  =  bJk  is 
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probability  of  classifying  successful  treatment  as  unsuccessful.  In  the  former  case 
the  task  is  discharged  correctly  as  completed;  in  the  latter  case  it  is  incorrectly 
re-served. 

•  b)k :  probability  of  successful/correct  assessment  as  unsuccessful  treatment  of  a 
task  of  type  j,  given  unsuccessfully  prosecuted/served  as  task  type  k.  In  this  case 
the  task  is  reclassified  and  served  again. 

•  There  are  two  customer  types  in  our  present  examples;  their  basic  service  times 
are  constant  (delta-function  distributed)  with  means  denoted  su  s2:  S1  =sx,  S2  =  s2 
with  probability  one. 

6.1  Discussion  of  the  Figures 

In  each  case  investigated,  we  display  the  probability  of  successfully  completing 
service:  reaching  the  server  through  the  initial  queue  and  subsequently  being  successfully 
served  before  deadline  elapse.  Three  models  are  compared:  J30  =  the  probability  of 
successful  completion  of  service  with  no  task  refusals  (5.3);  =  the  probability  of 

successful  completion  of  service  with  task  refusal  depending  on  the  virtual  waiting  time 
(D.ll);  and  /f  —  the  probability  of  successful  completion  of  service  with  task  refusal 
depending  on  the  virtual  waiting  time  and  allowed  service  time,  (DAS)  and  (D.ll).  We 
also  illustrate  the  numerical  quality  of  Approximation  II  in  the  graphs  that  follow,  and 
compare  Approximation  I  and  II  in  the  tables. 
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Prob.  of  successful  task  completion 


Probability  of  successful  task  completion 
Probability  of  correct  task  completion  assessment:  b=0.5  b*=0.5 
cii=0.7  ml  1=0.7  m22=0.7  0=0.1  s1=. 5  s2=1 


X:  Task  arrival  rate 

-• — Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time  and  allowed  service  time 

♦  -  -  Approximate  probability  of  successful  completion  of  service  with  refusal 
depending  on  the  virtual  waiting  time  and  allowed  service  time 

-a — Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time 

« -  -  Approximate  probability  of  successful  completion  of  service  with  refusal 
depending  on  the  virtual  waiting  time 

■* — Probability  of  successful  completion  of  service  for  system  with  no  refusal 

Figure  1 
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Table  for  Figure  1 


1 

task 

arrival 

rate 

prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

approx, 
prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

Approx  I/n 

prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

approx, 
prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

Approx  I/II 

prob.  of 
success, 
comple.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

approx, 
prob.  of 
success, 
comple.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

Approx  I/II 

prob.  of 
success, 
comple.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

approx, 
prob.  of 
success, 
compl.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

Approx  I/II 

prob.  of 
success, 
comple.  of 
service  for 
system 
with  no 
refusals 

0.2 

0.96 

0.96/0.96 

0.97 

0.97/0.97 

0.57 

0.57/0.57 

0.58 

0.58/0.58 

0.57 

0.91 

0.89/0.90 

0.92/0.92 

0.55 

CD 

O 

0.83 

0.79/0.81 

0.86 

0.84/0.85 

0.47/0.48 

0.51 

0.50/0.51 

0.31 

0.74 

0.68/0.70 

KBifli 

0.74/0.76 

0.44 

0.41/0.42 

0.47 

0.44/0.45 

0.00 

0.64 

0.58/0.60 

w§m 

0.64/0.66 

0.38 

0.35/0.36 

0.41 

0.38/0.39 

0.00 

1.2 

0.50/0.52 

Mflalii 

0.56/0.58 

0.33 

0.30/0.31 

0.36 

0.33/0.34 

0.00 

1.4 

0.48 

0.44/0.46 

0.53 

0.49/0.51 

0.28 

0.26/0.27 

0.32 

0.29/0.30 

0.00 

1.6 

0.42 

0.39/0.40 

0.47 

0.43/0.45 

0.25 

0.23/0.24 

0.28 

0.26/0.27 

0.00  1 

1.8 

0.37 

0.35/0.36 

0.41 

0.39/0.40 

0.22 

0.21/0.22 

0.25 

0.23/0.24 

0.00 

2.0 

0.33 

0.32/0.33 

0.37 

0.35/0.36 

0.20 

0.19/0.20 

0.22 

0.21/0.22 

0.00 

Parameters  for  Figure  1 

arrival  rate  of  tasks,  X  =  0.2  (0.2)  2 
prob.  of  correct  task  class:  cu  =  c22  =  0.7 

prob.  of  correctly  classifying  a  complete  task  as  complete:  bu  =  bn  =  b2l  =  b22  =  0.5 
prob.  of  correctly  classifying  an  incomplete  task  as  incomplete:  bn  =  b\2  =  b\x  —  b"22  =  0.5 
prob.  complete  task  of  type  j  that  is  correctly  classified  as  type  j:  mu=  m22  =  0.7 
prob.  complete  task  of  type  j  that  is  incorrectly  classified  as  type  k:  ml2  =  m2l  =  0 
service  time  for  task  classified  as  type  1:5,  =  0.5 
service  time  for  task  classified  as  type  2:  s2  =  1 
prob.  an  arriving  task  is  of  type  i :  px  =p2  =  0.5 
mean  of  the  exponential  deadline:  =  10.00 
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Discussion  of  Figure  1 

This  demonstrates  the  anticipated  decrease  in  the  probability  of  receiving  service 
(transiting  queue  before  deadline  elapse),  and  the  probability  of  ultimate  correct  task 
completion  as  the  arrival  rate,  A,  increases.  Notice  that  Approximation  II,  (5.16), 
relatively  closely,  but  conservatively,  tracks  the  exact  solution  of  (D.ll),  bounding  the 
latter  from  below.  Approximation  I  does  nearly  as  well.  The  payoff  from  being  able  to 
recognize  deadline  elapse  in  service  (and  task  ejection)  is  evident,  but  the  dramatic  effect 
is  caused  by  the  existence  of  a  deadline-recognized  queue  admission  capability:  if 
deadlines  are  ignored  then  the  queue  quickly  saturates  and  the  success  probability 
plummets.  This  happens  despite  the  fact  that  the  mean  deadline,  10.0,  is  much  greater 
than  the  mean  ISOT  (“service  time”),  1.5,  in  this  example. 

Next  we  study  the  effect  of  varying  the  classification  parameters. 
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Probability  of  successful  task  completion 
Probability  of  correct  task  completion  assessment:  b=0.9  b*=0.9 
ml  1=0.7  m22=0.7  6=0.1  X=.6  s1=.5  s2=1 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

c:  Probability  of  correct  task  classification 


— • —  Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time  and  allowed  service  time 

-  Approximate  probability  of  successful  completion  of  service  with  refusal 
depending  on  the  virtual  waiting  time  and  allowed  service  time 

— ■ — Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time 

-  -  *  -  -  Approximate  probability  of  successful  completion  of  service  with  refusal 
depending  on  the  virtual  waiting  time 

— * —  Probability  of  successful  completion  of  service  for  system  with  no  refusal 


Figure  2 
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Probability  of  successful  task  completion 
Probability  of  correct  task  completion  assessment:  b=0.9  b*=0.3 
ml  1=0.7  m22=0.7  0=0.1  si =.5  s2=1 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

c:  Probability  of  correct  task  classification 


— • — Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time  and  allowed  service  time 

Approximate  probability  of  successful  completion  of  service  with  refusal 
depending  on  the  virtual  waiting  time  and  allowed  service  time 

— * —  Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time 

-  •  *  -  -  Approximate  probability  of  successful  completion  of  service  with  refusal 
depending  on  the  virtual  waiting  time 

— * —  Probability  of  successful  completion  of  service  for  system  with  no  refusal 


Figure  3 
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Prob.  of  successful  task  completion 


Probability  of  successful  task  completion 
Probability  of  correct  task  completion  assessment:  b=0.3  b*=0.9 
ml  1=0.7  m22=0.7  0=0.1  ?i=.6  s1=. 5  s2=1 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

c:  Probability  of  correct  task  classification 


— • — Probability  of  successful  completion  of  service  with  refusal  depending  on 
the  virtual  waiting  time  and  allowed  service  time 

-  -  ♦  -  -  Approximate  probability  of  successful  completion  of  service  with  refusal 

depending  on  the  virtual  waiting  time  and  allowed  service  time 

— » — Probability  of  successful  completion  of  service  with  refusal  depending  on 
the  virtual  waiting  time 

-  -  •  -  -  Approximate  probability  of  successful  completion  of  service  with  refusal 

depending  on  the  virtual  waiting  time 

— * —  Probability  of  successful  completion  of  service  for  system  with  no  refusal 


Figure  4 
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Parameters  for  Figures  2-4 

prob.  of  correct  task  class:  cu  =  c22  =  0.1  (0.1)  1 

prob.  of  correctly  classifying  a  complete  task  as  complete:  bu  =  bn  =  b2X  =  b22  -  0.9 
(Figures  2  and  3),  0.3  (Figure  4) 

prob.  of  correctly  classifying  an  incomplete  task  as  incomplete:  b*n  =  bf2  =  b'2X  =  b22  =  0.3 
(Figure  3),  0.9  (Figures  2  and  4) 

prob.  complete  task  of  type  j  that  is  correctly  classified  as  type  j:  mu  =  m22  =  0.7 

prob.  complete  task  of  type  j  that  is  incorrectly  classified  as  type  k:  mn  =  m2x  =  0 

service  time  for  task  classified  as  type  1 :  sx  =  0.5 

service  time  for  task  classified  as  type  2:  s2  =  1 

arrival  rate  of  tasks,  A  =  0.6 

prob.  an  arriving  task  is  of  type  i:  px  =p2  =  0.5 

mean  of  the  exponential  deadline:  (<^''  =  10 


Table  for  Figure  2 


prob.  of 
correct 
task  class. 

prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

approx, 
prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

Approx  I/II 

prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

approx, 
prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

Approx  I/II 

prob.  of 
success, 
comple.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

approx, 
prob.  of 
success, 
comple.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

Approx  I/II 

prob.  of 
success, 
comple.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

approx, 
prob.  of 
success, 
compl.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

Approx  I/II 

prob.  of 
success, 
comple.  of 
service  for 
system 
with  no 
refusals 

0.1 

0.35 

0.31/0.34 

0.50 

0.44/0.47 

0.10 

0.09/0.10 

0.15 

0.13/0.14 

0.00 

0.2 

0.48 

0.42/0.45 

0.60 

0.53/0.56 

0.22 

0.19/0.21 

0.28 

0.25/0.26 

0.00 

0.3 

0.58 

0.52/0.55 

0.68 

0.62/0.65 

0.30/0.32 

0.39 

0.35/0.37 

wmm 

0.67 

0.61/0.64 

0.75 

0.69/0.72 

0.40/0.42 

0.49 

0.45/0.47 

w 

0.5 

— — 

0.69/0.71 

0.80 

0.75/0.77 

■H 

0.49/0.51 

0.57 

0.54/0.55 

M| 

0.6 

■ 

0.75/0.77 

0.84 

0.80/0.82 

0.57/0.59 

0.63 

0.61/0.62 

0.7 

0.84 

0.81/0.82 

0.87 

0.85/0.86 

0.64/0.65 

0.69 

0.67/0.68 

KB 

0.8 

0.87 

0.85/0.86 

0.89 

0.88/0.89 

0.72 

0.70/0.71 

0.74 

0.72/0.73 

0.62 

0.9 

0.90 

0.88/0.89 

0.91 

0.90/0.91 

0.76 

0.74/0.75 

0.77 

0.76/0.77 

0.70 

1.0 

0.91 

0.90/0.91 

0.93 

0.92/0.92 

0.79 

0.78/0.79 

0.79 

0.80/0.80 

0.76 
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Table  for  Figure  3 


prob.  of 
correct 
task  class. 

| 

prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  ■ 

approx, 
prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

Approx  I/II 

prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

approx, 
prob.  of 
receiving 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

Approx  I/II 

prob.  of 
success, 
comple.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time 

approx, 
prob.  of 
success, 
comple.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
'time 

Approx  I/II 

prob.  of 
success, 
comple.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

approx, 
prob.  of 
success, 
compl.  of 
service 
with 
refusal 
depending 
on  virtual 
waiting 
time  and 
allowed 
service 
time 

Approx  I/II 

prob.  of 
success, 
comple.  of 
service  for 
system 
with  no 
refusals 

0.1 

■39 

0.91/0.92 

0.94 

0.93/0.93 

0.08 

0.08/0.08 

0.08 

0.08/0.08 

0.08 

0.2 

■SEMI 

0.92/0.92 

0.94 

0.93/0.94 

0.16 

0.16/0.16 

0.16 

0.16/0.16 

0.15 

0.3 

0.93/0.93 

0.94 

0.94/0.94 

0.23 

0.23/0.23 
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Discussion  of  Figures  2-4 

These  figures  illustrate  the  possible  dependence  of  successfully  joining  the  queue,  and 
successful  task  completion  on  the  value  of  cn  =  c22  =  c:  the  probability  of  correctly 
classifying  the  task  type  (for  fixed  arrival  rate,  A,  and  other  parameters).  Apparently  there 
is  strong  dependence:  advantage  accrues  to  systems  with  task  refusals.  Accepting  all  tasks 
can  send  the  system  into  saturation  for  small  values  of  cn  =  c22  =  c;  it  only  approaches  the 
task  refusal  systems  if  c  closely  approaches  unity  (current  task  type  classification  is 
nearly  perfect).  The  size  of  the  advantage  depends  on  the  ability  to  correctly  assess  task 
completion.  If  the  probability  of  assessing  an  incomplete  task  as  incomplete  is  small  (b"  = 
0.3),  then  many  tasks  are  thrown  out  of  service  before  they  are  complete.  Thus,  the  traffic 
intensity  is  less  than  1  for  the  system  with  no  refusal  and  the  probability  of  successful 
task  completion  is  about  the  same  as  for  a  system  with  task  refusal.  If  the  probability  of 
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assessing  a  complete  job  as  complete  is  small  ( b  =  0.3),  the  additional  (non-productive) 
service  on  already  complete  tasks,  can  saturate  the  system  with  no  refusals;  it  also 
increases  the  probability  a  task  will  be  refused  in  a  system  with  refusal.  Thus,  increasing 
the  probability  of  correct  task  classification  can  have  less  effect  if  b *  is  small. 

In  the  next  case  we  explore  the  effect  of  quality  of  post-service  assessment. 
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Probability  of  successful  task  completion 
Probability  of  correct  task  completion  assessment:  b*=0.9 
cii=0.7  ml  1=0.7  m22=0.7  0=0.1  X=.6  si  =.5  s2=1 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

b:  Probability  of  assessing  incomplete  task  as  incomplete 


— • —  Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time  and  allowed  service  time 

-  Approximate  probability  of  successful  completion  of  service  with  refusal 
depending  on  the  virtual  waiting  time  and  allowed  service  time 

— » — Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time 

-  -  *  -  -  Approximate  probability  of  successful  completion  of  service  with  refusal 
depending  on  the  virtual  waiting  time 

— * —  Probability  of  successful  completion  of  service  for  system  with  no  refusal 


Figure  5 


b.  of  successful 


-Probability  of  successful  completion  of  service  with  refusal  depending  on 
the  virtual  waiting  time 

-Approximate  probability  of  successful  completion  of  service  with  refusal 
depending  on  the  virtual  waiting  time 

-  Probability  of  successful  completion  of  service  for  system  with  no  refusal 
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Parameters  for  Figures  5-7 

prob.  of  correctly  assessing  a  complete  task  as  complete:  bn  =  bn  =  b21  =  b22  =  0.1  (0.1)  1 
(Figure  5),  bu  =  bl2  =  621  =  622  =  0.9  (Figure  6),  bu  =  b12  =  b2l  =b22  =  0.1  (Figure  7) 
prob.  of  correct  task  class:  cu  =  c22  =  0.7 

prob.  of  correctly  assessing  an  incomplete  task  as  incomplete:  bu  =  b*2  =  b2]  =  b22  =  0.9 
(Figure  5),  bu  =  b,2  =  b21  =  b22  =  0.1  (0.1)  1  (Figure  6),  b^  =  b^2  =  b21  =  b22  =  0.7 
(Figure  7) 

prob.  complete  task  of  type  j  that  is  correctly  classified  as  type  j:  mu  =  m22  =  0.7 

prob.  complete  task  of  type  j  that  is  incorrectly  classified  as  type  k:  mn  =  m2l  =  0 

service  time  for  task  classified  as  type  1 :  s,  =  0.5 

service  time  for  task  classified  as  type  2:  s2  =  1 

arrival  rate  of  tasks,  A.  =  0.6  (Figures  5-6),  A  —  0.8  (Figure  7) 

prob.  an  arriving  task  is  of  type  i:  p{  =p2  =  0.5 

mean  of  the  exponential  deadline:  (0)'*  =  10  (Figures  5-6),  0=  (0.1  (0.2)  (1.9))  (Figure  7) 

Discussion  of  Figures  5-6 

Increasing  the  probability  of  correctly  assessing  an  incomplete  task  as  incomplete 
by  =  b"  results  in  increases  in  the  probability  of  correct  task  completion  for  those  systems 
with  refusals.  However,  for  the  system,  with  no  refusals,  increasing  by  =  b‘  results  in 
larger  service  times  and  thus  decreases  the  probability  of  successful  task  completion. 

Increasing  the  probability  of  correctly  assessing  a  complete  task  as  complete,  by  =  b 
results  in  increases  in  the  probability  of  correct  task  completion  for  those  systems  with 
and  without  customer  refusals.  For  b  <  0.7,  the  system  with  no  refusals  is  saturated  and 
the  probability  of  correct  task  completion  is  0. 
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Probability  of  successful  task  completion 
Probability  of  correct  task  completion  assessment:  b=0.1  b*=0.7 
cii=0.7  ml  1=0.7  m22=0.7  ^=0.8  si  =5  s2=1 


— • —  Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time  and  allowed  service  time 

-  -  ♦  -  -  Approximate  probability  of  successful  completion  of  service  with  refusal 

depending  on  the  virtual  waiting  time  and  allowed  service  time 

— « — Probability  of  successful  completion  of  service  with  refusal  depending  on  the 
virtual  waiting  time 

-  - » -  •  Approximate  probability  of  successful  completion  of  service  with  refusal 

depending  on  the  virtual  waiting  time 

— * —  Probability  of  successful  completion  of  service  for  system  with  no  refusal 


Figure  7 
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Discussion  of  Figure  7 

Figure  7  displays  the  probability  of  successful  service  completion  as  a  function  of  the 
deadline  rate,  6.  In  Figure  7  the  probability  of  correctly  assessing  a  complete  job  as 
complete  is  small;  b  =  by-  0.1.  Note  that  for  the  case  of  task  refusal  depending  on  both 
the  virtual  waiting  time  and  the  allowed  time  in  service,  there  is  an  optimum  6.  Reason: 
the  small  probability  of  assessing  a  complete  task  as  being  complete  is  creating  much 
unproductive  work  for  the  server.  A  <9  that  is  too  small  can  result  in  a  task  remaining  in 
service  when  it  is  complete.  A  9  that  is  too  large  results  in  too  many  tasks  being  turned 
away. 

7.  Summary 

A  stochastic  model  has  been  introduced  that  allows  initial  discussion  of  the  general 
problem  of  uncertain  impatient  service.  The  influence  of  the  various  processes  that  may 
affect  such  service  has  been  numerically  explored  in  special  cases;  it  can  be  seen  that  the 
ability  to  adapt  to  long  queues  (by  refusing  admission)  as  in  Section  5.2  and  Appendix  D, 
can  improve  overall  performance.  In  fact,  if  deadlines  did  not  exist  they  might  well  be 
imposed  in  order  to  improve  overall  long-run  system  performance.  Likewise, 
improvement  of  performance  by  improving  either  pre-service  classification,  and/or  post¬ 
service  assessment  can  be  quantitatively  traced  and  choices  made.  Degradation  of  either 
of  the  latter  capabilities  can  substantially  degrade  overall  system  performance,  as 
measured  by  the  probability  of  successful  complete  task  processing  (before  deadline 
elapse). 

The  present  paper  scratches  the  surface  of  an  important  and  largely  neglected  service 
system  design  and  control  problem.  It  is  planned  to  pursue  other  ramifications,  such  as 
non-stationary  phenomena  (via  fluid  approximations)  and  adaptive  control  (dynamic 
priorities)  in  subsequent  work. 
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APPENDIX  A 
Second  Moment  of  ISOT 

Squaring  and  taking  conditional  expectations  leads  to  an  expression  for  the  second 
moment  of  completion  time  (needed  for  calculating  long-run  expected  system  occupancy 
when  tasks  are  queued).  We  express  the  formula  in  terms  of  the  expressions  for  the  mean 
E[Cj\,  (4.4),  and  the  mean  and  second  moment  of  recognition  time  Kf. 

£[c;]  =  |X^£[S1!]+22;cJ»£[Sym}.(S,)>;£[C,] 

k 

k  k 

1  ^  cjkE\mjk  (Sk 

k 

where 

(a.2) 

and 

4jr']=(Sc/»£[s*]+2Ec^‘]^£[Ji)]]-i-[1-ZcAj-  (a.3) 


(A.1) 
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APPENDIX  B 

Model  with  Unobservable  Exponential  Deadlines 


The  Laplace-Stieltjes  transform  (LST)  of  the  “completion”  time  can  be  calculated  by 
taking  conditional  expectations  in  (4.1).  The  result  is 

]  =  {s  )£>  ]] + X  ^  ]} 


*■  1  -  X  mjk  (Sk  )\b]k 


E\e'Xj  ]  =  \  J]  C;*£[<T®‘  J&yit  [  +  1  -  2  ]b 


Note  that  the  above  can  be  interpreted  as  the  probability  that  a  sojourn  at  the  server, 
including  repeats,  ends  before  the  termination  of  an  independent  exponential  random 
variable  with  mean  1  Id. 

To  obtain  the  probability  that  a  task  service  has  been  satisfactorily  completed  at 

sojourn  completion  and  that  sojourn  ends  before  the  termination  of  the  independent 

exponential  random  variable  with  mean  1 19,  define 

Sk  with  probability  cjkmjk  (Sk ); 

dj=\  _  (B.3) 

Sk  +Dj  with  probability  cjkmjk(Sk)bjk. 

This  is  simply  (4.1)  with  a  term  omitted.  It  is  seen  that 


E\e  ®']  =  k  1  - £ mjk (Sk )b]k ] 
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APPENDIX  C 

A  Model  with  the  Option  to  Reclassify 
Unobservable  Deadlines 


In  this  appendix  we  present  results  for  a  model  in  which  a  server  may  opt  to  reclassify 
the  task  type  if  it  perceives  that  the  task  is  not  complete.  There  is  a  time  penalty  A  for 
reclassification.  Let  or  be  the  probability  that  the  server  decides  to  reclassify  the  task  after 
a  task  service  which  is  perceived  not  to  have  completed  the  task. 


'sk 

with  probability  [mjk  (Sk  )bjk  +  mjk  (Sk  )bjk  ]; 

Sk  +  Cjk 

with  probability  mJk  ( Sk  )b]k  (1  -  a); 

Sk  +  Cj  +  A 

with  probability  m )k  (Sk )b*k a; 

Sk  4"  Kjk 

with  probability  m  jk  ( Sk  )bjk  (1  -  a); 

Sk  +  Kj 

with  probability  mjk  ( Sk  )bJk  a\ 

CJk  represents  the  service  time  of  a  task  of  type  j  that  has  been  classified  as  a  type  k.  Cj 
represents  the  service  time  of  a  task  of  type  j.  Kjk  represents  the  random  time  until  a 
completed  task  of  type  j  that  has  been  classified  as  a  type  k  is  so  identified.  X}  represents 
the  random  time  until  a  completed  task  of  type  j  is  so  identified,  the  recognition  time: 

f Sk  with  probability  b Jk 


Kjk 


•  Sk  +  K'Jk  with  probability  bjk  (1  -  a) 
Sk  +  Kj  +  A  with  probability  bjka 


(C.2) 


Cj  =  Cjk  with  probability  Cjk ,  k  =  1, . . . ,  J 
Kj  =  Kjk with  probability  c jk ,  k  =  1, . . . ,  J 


Notice  that  the  task  accomplishment  probability  is  allowed  to  depend  explicitly  on  the 
allocated  service  time,  Sk,  and  that  the  “completion”  time  can  terminate  with 
unrecognized  incomplete  task  service,  i.e.  incompletely.  The  random  variables  C'-  ,C'Jk 
JK) ,  and  K'Jk  above  are  independent  stochastic  replicas  of  Cj  CJk  A)  and  Kjk,  as  usual. 
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The  expectations  required  are  now  calculated. 


Thus 


where 


Similiarly 


where 


Further, 


and 


1 ]*>,*« 


k 


AW- 


£[e’®]6 


7* 


l-£[e-ft]6j*(l-ar) 


E[e-mt]E{e-*]bjka 
~  \-E[e-^]bjk(l- a) 


ATyt(C)+&(04-*'] 

[  1_  DAC) 

Njic  (C)  =  [m  jk  (Sk  )bJk  +  mjk  (Sk  )bjk  ]] 

+E[e'^  ]E[e'*  ]£[e-*'  ]bJka 
+E[e~mk  1 njk  ( Sk  )^e'Xjk  ] mjkbjk  (1  -  a). 

Djk{C)  =  1  -  E[e-^mjk{Sk)f;a 


(C.3) 


(C.4) 


(C.5) 


(C.6) 


(C.7) 


(C.8) 
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After  appropriate  weighting  by  pp  the  probability  that  an  arrival  is  of  type  j,  the  needed 
moments  and  transforms  can  be  calculated  as  was  done  previously. 
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APPENDIX  D 

Forward  Kolmogorov  (Taka^s-Benes)  Equations  with  Exponential  Balking 

Suppose  that  tasks  arrive  at  a  service  facility  according  to  a  Poisson  process  with  rate 
A.  Service  times  are  independent  and  identically  distributed.  Let  W(t)  be  the  total  virtual 
work  in  the  system  at  time  t.  Each  task  has  a  deadline  which  is  exponentially  distributed 
with  mean  1/6  if  the  waiting  time  or  virtual  work  present  when  the  task  arrives  exceeds 
the  deadline  the  task  does  not  enter  the  system.  This  approximates  the  situation  in  which 
tasks  whose  deadlines  have  elapsed  when  they  reach  the  server  are  not  served.  With  some 
modification  it  addresses  the  situation  in  which  a  deadline  elapses  during  service. 

D.l  Statistically  Identified  Deadlines  and  Service  to  Completion 

Let  the  distribution  function  of  W(t)  be 

Fw{x;t;0)  =  P{W{t)<x} 

and  express  this  as 

x 

Fw{x\  t;  0 )  =  p0  ( t ;  0) + J  p(z;t;  0)dz, 

0 

where 


po{t;0)  =  P{W{t)  =  O}. 

Since,  given  W(t),  the  task  joins  the  queue  with  probability  e~^‘\  the  probability  its 
deadline  does  not  expire  while  in  queue,  one  can  write 


p(x;t + A t;0)  =  p(x  +  A  t;t;  6)[l  -  Ae'{x+6l)ffAtj 

+ p0  (V;  0)b(x)AAt  +  AAt  [  e~®p(y;  t;  0)b(x  -  y)dy  +  o(At) 

JQ 

where  b  is  the  density  function  of  the  positive  service  time  C.  Also, 


(D.la) 


Po  (t  +  At;  0)  =  p0  (t;  6)[l  -  AAt] + p(  0,  t;  0)At(l  -  AAt)  +  o(At) .  (D.  1  b) 
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Taking  limits  as  At  -»  0  results  in 

4~Po  &  0)  =  ~M  {*', 9)  +  P(0,  t;  6) , 

at 

and  also 

rl 

— p(x;  t;0)  =  —  p(x;  t;  9)  -  p(x;  t ;  0)Xe~6c  +  Xp0  ( t ;  0)b{x) 
ot  ox 

+Xfe~^p{y;t;0)b(x -y)dy 

JO 

If  t  -» oo,  then  a  steady-state  density  satisfies 

0  =  —p(x;  9)-p(x;  0)Ae~&  +  p0b(x)X  +  X  [  e~® p(y)b(x  -  y)dy , 
ox  J0 

0  =  -Xp0  +  p{  0). 

Laplace  transform  to  obtain 

p  (s;  0)  =  Jo+  e~sxp(x;  0)dx  and  b*  (s)  =  Jq  e~sxb(x)dx . 

Thus,  (D.4a)  implies 


Let 


then 


where 


sp"  (5;  0)  =  [p{  0;  0)  +  Xp‘  (s  +  0, 0)]  [l  -  b*  (5)] 

p{s;  0)  =  f  e~sxdFw(x;  0)  =  p0  ( 9 )  +  JT  esxp{x\  0)dx 
=  po(0)+p\s\0)-, 

y/{s\  0)  =  po{0)+ py/{s  +  0, 0)S{s) 


0{s)  = 


l-b*{s) 
sE[C]  ' 


Since  i/X 0;  9)  =  1  =p0  +  pyX.9  9), 


(D.2) 


(P-3) 


(D.4a) 


(D.4b) 


(D.5) 


(P-6) 


(P-7) 


(P-8) 
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po=l-pv{&,0)  =  \-M[e-av]E[C\.  (D.9) 

which  motivates  the  heuristic  approximation  (actually  lower  bound)  of  Section  5.2. 

Iterative  solution  to  the  equation  (D.7) 

Since 

0)  =  [l -py/(&,  0)] + pyAs  +  0 ;  0)S(s) 
it  follows  that,  putting  s  =  0,  and  defining 

y/{0)  =  y/(0,  0)  =  [\-  pip(0)] + py(20, 0)0(0) 
and 


ip{20)  =  y/(20, 0)  =  [l  -  pip{0)\ + pys($0,  0)S(20) . 
Substituting  the  expression  for  y{20)  into  that  for  y^0)  results  in 

y(0)  =  [l-  pyA0)] +/?{[!-  P>P(0)] + pip(20-,  0)0(20, 0)}S(0) 
=  l+p6(0)-pip(0)[l+p6(0)]+/0ip(30,0)6(20)6(0). 
Continuing  in  this  maruier  results  in  the  equation 

<p(0)  =  ^4(0, 0) + . . . + A(n;  0)  -  py/(0^A(O, 0 )+...  +  A(n ;  0)\ 
+py/((n  +  \)0, 0)A(n;  0). 

where 


(D.lOa) 


(D.lOb) 


A(0, 0)  =  1 

A(n ;  0)  =  pnS(n  0)S((n  - 1)<?)  x . . .  x  6(0). 

For  0>  0,  A(n;0)  ->•  0.  Thus  the  probability  that  an  arriving  task  joins  and  survives  the 
queue  before  deadline  elapse  is 


'£A(k;ff) 

- ■  (DU) 

1  +flYJA(k-,0) 

*=0 
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It  is  clear  that  the  infinite  sums  converge  exponentially  rapidly  for  0>  0,  and  that  this  is 
true  for  any  /> value. 

Further,  for  s  *  0 

*Ki;<?)  =  [l P-12) 

k= 0 

where 

C(k;s)  =  /  fj  S{s  +  W),  k  >  1  (D.13) 

»=o 

and 


C(0;s)=l. 


D.2  Services  Subject  to  Exponential  Deadline 

If  a  task  deadline’s  elapse  is  detectable  during  service  and  the  task  then  removed,  then 
the  distribution  of  service  time,  C,  must  be  replaced  by  that  of  CT=  min(C,  deadline),  the 
allowed  service  time.  Consequently  the  service  times  that  contribute  to  the  virtual 
waiting  time  are,  thanks  to  the  exponential  deadline  assumption,  iid  with  mean 


and  tail-transform  now 


(D.14) 


(0+s)£[Ct]  S(0) 

These  replace  E[C\  in  p,  and  in  the  previous  solution,  (D.l  1). 


(D.15) 


D.3  Class-Specific  Deadlines 

A  natural  generalization  of  the  above  is  to  allow  independent  Poisson  arrivals  from  J 
task  classes,  the  f'  rate  being  Aj  with  service  time  density  bj  and  exponential  deadline 
parameter  0y  Arguments  analogous  to  those  in  Appendix  D  enable  us  to  show  that 
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(D- 16) 


(D-17) 


(D.18) 


The  equation  (D.16)  can  be  solved  in  closed  form  (a  series)  by  successive 
substitution/iteration,  but  this  step  is  omitted. 
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